Multifractality of self-avoiding walks on percolation clusters 
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We consider self-avoiding walks (SAWs) on the backbone of percolation clusters in space dimen- 
sions d = 2, 3, 4. Applying numerical simulations, we show that the whole multifractal spectrum of 
singularities emerges in exploring the peculiarities of the model. We obtain estimates for the set of 
critical exponents, that govern scaling laws of higher moments of the distribution of percolation clus- 
ter sites visited by SAWs, in a good correspondence with an appropriately summed field-theoretical 
£ = 6 - d-expansion (H.-K. Janssen and O. StenuU, Phys. Rev. E 75, 020801(R) (2007)). 
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When studying physical processes on complicated frac- 
tal objects, one often encounters the interesting situation 
of coexistence of a family of singularities, each associated 
with a set of different fractal dimensions 1]. In these 
problems, the conventional scaling approach cannot de- 
scribe the system. Instead, an infinite set of critical ex- 
ponents is needed to characterize the different moments 
of the distribution of observables, which scale indepen- 
dently. These peculiarities are usually referred to as mul- 
tifractality [4I . The multifractal spectrum can be used to 
provide information on the subtle geometrical properties 
of a fractal object, which cannot be fully described by 
its fractal dimensionality. Indeed, clusters generated by 
diffusion-limited aggregation (DLA) [3] and percolation 
clusters have the same fractal dimensions, but a com- 
pletely different geometrical structure, which can be clar- 
ified, e.g., by studying the multifractality of the voltage 
distribution in percolation clusters j^] and the growth 
probability distribution in DLA [5|, |g] . Multifractal prop- 
erties arise also in many different contexts, for example in 
studies of turbulence in chaotic dynamical systems and 
strange attractors [2, LZ| , human heartbeat dynamics [8| , 
Anderson localization transition [9] etc. 

To understand the common roots underlying this phe- 
nomenon, it is worthwhile to consider the generic case of 
self-avoiding walks (SAWs) on fractal clusters. It is well 
established that configurational properties of SAWs on a 
regular lattice are governed by scaling laws; e.g. for the 
averaged end-to-end distance (r) of a SAW with N steps 
one finds in the asymptotic limit N —^ 00: 



(r) - N" 



(1) 



where i^sAW is an universal exponent depending only on 
the space dimension d. 

The scaling of SAWs changes crucially, when the un- 
derlying lattice has a fractal structure. Indeed, new crit- 
ical exponents were found for SAWs residing, e.g., on a 
Sierpinski gasket and Sierpinski carpet [10|. A related 
problem arises when studying SAWs on disordered lat- 



TABLE I: Fractal dimensions of percolation cluster dp^ and 
backbone of the percolation cluster dp^ , and correlation length 
critical exponent Up for different space dimensions d. 



d 



dt. 



df. 



2 91/49 [21J 1.650 ± 0.005 [23] 4/3 [24] 

3 2.51 ± 0.02 [22] 1.86 ± 0.01 [23] 0.875 ± 0.008 [25] 

4 3.05 ± 0.05 [22] 1.95 ± 0.05 [23] 0.69 ± 0.05 [23] 



tices with concentration p of structural defects very close 
to the percolation threshold pc- In this case, an incipient 
cluster of pure sites can be found in the system. The 
diameter of a typical cluster below Pc is characterized by 
the correlation length ^, which diverges as ^ ^ {p—Pc)~'^'' 
with an universal exponent Vp. Note that percolation 
clusters are fractal objects (see Table 1) and apparently 
change the universality class of residing SAWs; the scal- 
ing (1) holds in this case with an exponent Vp^ ^ t'sAW- 
Aiming to study the scaling of SAWs on a percolative lat- 
tice, we are interested rather in the backbone of percola- 
tion cluster: the structure left when all "dangling ends" 
are eliminated from the cluster. Infinitely long chains can 
only exist on the backbone of the cluster. 

Although the behavior of SAWs on percolative lattices 
served as a subject of numerous numerical [ijj, [l^, llS], 
3, [H [Hi and analytical [13, El, \lk, M] studies since 
the early 80th, not enough attention has been paid to 
clarifying the multifractality of the problem. Following 
an early idea of Meir and Harris 11211. it was only recently 
proven in field-theoretical studies [19|, |20| that the expo- 



nent Vp^ alone is not sufficient to completely describe the 
peculiarities of SAWs on percolation clusters. Instead, 
a whole spectrum ly^'^' of multifractal exponents emerges 
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with e = 6 — d. Note that putting q = in ^, we 
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FIG. 1: Different SAW trajectories with fixed end-to-end dis- 
tance R on the backbone of a percolation cluster in d = 2. 



restore an estimate for the dimension dp^ of the underly- 
ing backbone of percolation clusters found previously (26| 
yia J^*^"' = 1/dp^, whereas q ^ oo gives the percolation 
correlation length exponent zy(°°) = i/p^ which restores, 
e.g., the result of Amit [27|. f^*-^-* gives us the exponent 
i/p^ , governing the scaling law for the averaged end-to-end 
distance of SAWs on the backbone of percolation clus- 
ters. This is a major step forward, but the reliability of a 
two-loop e-expansion in e = 6 — d remains a priori ques- 
tionable in physical dimensions d — 2,3 where £ is large. 
To settle this question, we report in this letter a care- 
ful computer simulation study of SAWs on percolation 
clusters. We primarily aim to obtain precise numerical 
estimates for the multifractal exponents i'^'^' in a wide 
range of q and find the related spectrum of singularities 
on the underlying fractal cluster, thus going into deeper 
understanding of peculiarities of the model. 

We consider site percolation on regular lattices of edge 
lengths up to Lmax=400, 200, 50 in dimensions d—2, 3, 4, 
respectively. Each site of the lattice is occupied ran- 
domly with probability pc and empty otherwise. If for 
a given lattice it is not possible to find a spanning clus- 
ter, this disordered lattice is rejected and a new one is 
constructed. To extract the backbone of a percolation 
cluster, we apply the algorithm proposed in Ref. 28] . 
First, the starting point - "seed" - is chosen at the cen- 
ter of the cluster. Then, for all the sites on the borders 
of the lattice, the shortest paths to the "seed" are found, 
forming the so-called elastic backbone [29|. Finally, con- 
sidering successively each site of the elastic backbone and 
checking whether a "loop" (path of sites connected to the 
elastic backbone in two places at least) starts on it, we 
obtain the geometric backbone of the cluster. We con- 
struct 1000 clusters in each space dimension. 

Starting on the "seed" of a single cluster, we construct 
a SAW on it, applying the pruned-enriched chain-growth 
algorithm |30]. We let a trajectory of SAW grow step by 
step, until it reaches some prescribed distance (say R) 
from the starting point. Then, the algorithm is stopped, 
and a new SAW grows from the same starting point. In 



FIG. 2: Number of sites n{'w) of the backbone of a percolation 
cluster with weights w, visited by SAWs with fixed end-to-end 
distance R = 160, in d = 2 dimension. 



such a way, we are interested in constructing different 
possible trajectories with fixed end-to-end distance, as is 
shown schematically in Fig. [TJ For each lattice size L, we 
change i? up to « i/3 to avoid finite-size effects, since 
close to lattice borders the structure of the backbone of 
percolation clusters is distorted and thus can falsify the 
SAW statistics. 

Let us denote by K{R) the total number of constructed 
SAW trajectories between and R (we perform ~ 10^ 
SAWs for each value of R). Then, for each site i of the 
backbone we sum up the portion of trajectories, pass- 
ing through this site. In such a way, we prescribe a 
weight w{i) — K{i)/K{R) to each site i G R, and thus re- 
ceive some "population" , occupying the underlying frac- 
tal cluster. The distribution n{'w) of weights, prescribed 
to the sites of the backbone of a percolation cluster, vis- 
ited by SAWs with fixed end-to-end distance R — 160 in 
two dimensions, is shown in Fig. [2l 

The multifractal moments M^"^^ are defined as follows: 



M^^) 



E^'W 



(3) 



i£R 



Averaged over different configurations of the constructed 
backbones of percolation clusters, they scale as: 



Mil) ^ R^/'^'- 



(4) 



with exponents i/i'^' that do not depend on q in a linear or 
affine fashion, implying that SAWs on percolation clus- 
ters are multifractals. To estimate the numerical values 
of I'i'^) on the basis of data obtained by us (see Fig. [3]), 
linear least-square fits are used. The x^ value (sum of 
squares of normalized deviation from the regression line) 
serves as a test of the goodness of fit. 

It is clear that at q — we just count the num- 
ber of sites of the cluster of linear size R, and thus 
l/iyi^) corresponds to the fractal dimension of the back- 
bone d^^. Our results give d^^{d=2) = 1.647 ± 0.006, 
d^^{d=Z) = 1.865 ± 0.006, d^^{d=A) = 1.946 ± 0.006, in 
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FIG. 3: Averaged moments Mt?) as function of R in dou- 
ble logarithmic scale in d = 2, 3, 4. In each case, q — 0, 
1, 2, 3, 4, . . . going from above. Lines are guides to the eyes. 



very good agreement with Table 1. At q = 1, we restore 
the value of the exponent Vp^ , governing the scaling law of 
the end-to-end distance for S AWs on the backbone of per- 
colation clusters. We obtain iy'^^\d=2) = 0.779 ± 0.006, 
lyW (d=3) = 0.669 ± 0.006, ly^^'' (d=4) = 0.591 ± 0.006, in 
perfect agreement with our recent estimates [la] based 
on the scaling of the end-to-end distance with the num- 
ber of SAW steps. At g ^ oo, the so-called "red sites" of 
the backbone are mainly taken into account - the singly 
connected sites, such that cutting off one of them will 
produce a disconnection of the cluster. These sites are 
most often visited by SAWs, and thus have the maxi- 
mum weights. It was proven 3l|, that the number of 
"red sites" scales with linear distance i? of a cluster as: 
^rcd(-R) ^ R^/^^. Indeed, as it follows from our data, 
at large q the value of the exponent i^'''' tends to Vp^ 
the percolation correlation length exponent, cf. Table [H 
So, in two limiting cases (5 — * 0, g — > 00) we restore 
the corresponding multifractal exponents of the voltage 
distribution on the backbone of percolation clusters [4| . 
Note that in the latter problem, the exponent governing 
the scaling of the moment with q= 1 gives the resistance 



(1)^ 



exponent: e.g., v)^' {d=2) ~ 0.97 



FIG. 4: Left: spectrum of multifractal exponents v^'^' as func- 
tion of g in d = 2,3,4, dotted lines present [l]/[2] Fade ap- 
proximants to the analytical results of Janssen and StenuU 
(Eq. ^). Right: Spectral function /(a) in d = 2,3,4, the 
maximum value of /(a) gives the fractal dimension of the 
underlying backbone of percolation cluster. 



Due to the long tail of the distribution n{w), the pre- 
cision of our estimates decreases with increasing q. This 
problem turns out to be especially crucial when explor- 
ing the moments with negative powers q: the sites with 
small probabilities to be visited, which are determinant 
in negative moments, are very difficult to probe. 

Our estimates of the exponents i'^'^' for different q are 
presented in the left panel of Fig. [H These values appear 
to be in an astonishingly perfect correspondence with an- 
alytical estimates down to d = 2 dimensions, derived by 
applying Fade approximation to the e = 6 — d-expansion 
(P|, presenting the given series as ratio [r7T,]/[r7,] of two 
polynomials of degree m and n in e. We used the [l]/[2] 
approximant, because it appears to be most reliable in 
restoring the known estimates in limiting cases {q = 0, 
q — > 00). A direct use of the expression ([2]) gives worse 
results, especially for low dimensions d where the expan- 
sion parameter e = 6 — d is large. 

It is well known that the set of exponents governing 
scaling of multifractal moments of the type (JH) is related 



to the spectrum of singularities /(a) of fractal measure 
[5|, called also the spectral function. The physical mean- 
ing of /(a) in our problem is that the number Nfi(a) of 
sites i, where the weight w(i) scales as R~", behaves as: 



Nnia) r^ Rf'-'^\ 



(5) 



The singularity spectrum /(a) is given by the Legendre 
transform: 



f{a) ^ qa-T{q), a{q) 



dr(g) 
dq ■ 



(6) 



with T{q) = l/v'^'^\ Spectral functions f[a), obtained 
on the basis of our results, are given in the right panel 
of Fig. m The general properties of /(a) are as fol- 
lows: it is positive on an interval [ckmin , ctmax] , where 



= lim T{q)/{q 



!)■ 



amin = lim T{q)/{q - 1), ar,_____ 

q — *+oc q — ^ — u*j 

The maximum value of the spectral function gives the 
fractal dimension of the underlying structure, which in 
our case corresponds to the dimension of the backbone 
of percolation clusters. 

To conclude, we have shown numerically that SAWs 
residing on the backbone of percolation clusters give rise 
to a whole spectrum of singularities, thus revealing multi- 
fractal properties. To completely describe peculiarities of 
the model, the multifractal scaling should be taken into 
account. We have found estimates for the exponents, 
governing different moments of the weight distribution, 
which scale independently, in surprisingly good coinci- 
dence with two-loop £-expansions. The behaviour of the 
spectral function, describing the frequency of observa- 
tion of a set of singularities on the underlying backbone 
of percolation clusters, is analyzed as well. 
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